Relaxation Dynamics and Interrupted Coarsening 
in Irrationally Frustrated Superconducting Arrays 
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Equilibrium and non-equilibrium relaxation behaviors of two-dimensional superconducting arrays 
are investigated via numerical simulations at low temperatures in the presence of incommensurate 
transverse magnetic fields, with frustration parameter / = (3 — Vo)/2. We find that the non- 
equilibrium relaxation, beginning with random initial states quenched to low temperatures, exhibits 
a three-stage relaxation of chirality autocorrelations. At the early stage, the relaxation is found to 
be described by the von Schweidler form. Then it exhibits power-law behavior in the intermediate 
time scale and faster decay in the long-time limit, which together can be fitted to the Ogielski form; 
for longer waiting times, this crosses over to a stretched exponential form. We argue that the power- 
law behavior in the intermediate time scale may be understood as a consequence of the coarsening 
behavior, leading to the local vortex order corresponding to / = 2/5 ground-state configurations. 
High mobility of the vortices in the domain boundaries, generating slow wandering motion of the 
domain walls, may provide mechanism of dynamic heterogeneity and account for the long-time 
stretched exponential relaxation behavior. It is expected that such meandering fluctuations of the 
low-temperature structure give rise to finite resistivity at those low temperatures; this appears 
consistent with the zero-temperature resistive transition in the limit of irrational frustration. 

PACS numbers: 74.50+r, 67.40.Fd 



I. INTRODUCTION 

Relaxation properties of systems with a great num- 
ber of metastable states have attracted much atten- 
tion in the recent decades These systems arc 
usually characterized by the existence of both disorder 
and frustration. A well-known system is the Ising spin 
glass modelji which displays non-exponential relaxation 
at low temperatures^^ Interesting aging phenomena 4 
were also observed in this system^ The Coulomb glass 
model, which is another interesting disordered system, 
has recently been shown to exhibit relaxation of the 
stretched exponential form. 6 Such non-exponential relax- 
ation was also reported in frustrated systems without dis- 
order and its relation with the percolation transition was 
discussed* 7 - 

In a superconducting array, frustration can be induced 
in a controllable way by applying an external magnetic 
field. It has crucial effects on thermodynamics of the sys- 
tem and results in a variety of equilibrium properties. 8,9 
Such remarkable diversity in the equilibrium properties 
naturally leads to expectations that a variety of inter- 
esting dynamic behaviors is also present, particularly in 
relaxation toward equilibrium. 

An interesting limiting situation arises in the presence 
of irrational frustration, the most typical case of which 
is provided by the frustration parameter / = 1 — g with 
the golden number g = (s/E— l)/2. The system was first 
suggested to exhibit a spin-glass-like phase at low tem- 
peratures due to self-generated disorder.— Subsequently, 
it was argued to display novel finite-size effects^ 1 - that the 
size-dependent transition temperature decreases mono- 



tonically with the system size, resulting in the absence 
of a finite-temperature transition^ Simulation results 
for the current- voltage characteristics supported the zero- 
temperature transition^ while experimental results were 
interpreted to exhibit finite-temperature transitions! 14 ! 15 
Recent Monte Carlo simulations, investigating the vortex 
configuration at low temperatures in the system with / 
given by rational approximants to 1— g, indicate the exis- 
tence of a low-temperature phase where the helicity mod- 
ulus takes a finite value along one direction and vanishes 
along the other down to very low temperatures . 16 ' 17 ! 18 
From dynamical perspectives, due to the existence of 
many metastable states that are almost degenerate with 
one another (which is attributable to the incommen- 
surate magnetic field), one may expect characteristic 
slow relaxation in the array with irrational frustration. 
Dynamic simulations of this system, based on simple 
Langevin dynamics, indeed disclosed a crossover tem- 
perature below which strongly nonexponential relaxation 
emerges, exhibiting some analogy to the behavior of su- 
percooled liquids^ 

In this paper we investigate the relaxation behavior 
of the superconducting array with irrational frustration 
in both equilibrium and nonequilibrium situations, em- 
ploying the resistively-shunted-junction (RSJ) dynamics 
in the overdamped limit (i.e., junction capcitances are 
neglected). Note that the present RSJ dynamics can, 
in principle, be realized in real junction array experi- 
ments; this is in contrast to the Langevin dynamics em- 
ployed in existing simulations, which assumes hypothetic 
dissipation between the superconducting islands and the 
ground. 
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At equilibrium, we observe that chirality autocorrela- 
tions are characterized by stretched exponential relax- 
ation (with a temperature-dependent stretching expo- 
nent) in a wide range of intermediate and low temper- 
atures. The relaxation time exhibits non-Arrhenius be- 
havior with the Vogel-Tammann-Fulcher type of temper- 
ature dependence. 

On the other hand, beginning with random initial 
states quenched to low temperatures (T < 0.14), the dy- 
namics exhibits slow aging behavior, not reaching the 
equilibrium relaxation within our computing time. It 
is observed that, for a short time (approximately up to 
time t w 10 to 20 depending on the waiting time and 
temperature), relaxation of the chirality autocorrelation 
function exhibits the von Schweidler behavior^ After 
the short time, chirality autocorrelations for short wait- 
ing times relax according to the so-called Ogielski form 
with temperature-dependent exponents. For longer wait- 
ing times, this behavior slowly changes into the stretched 
exponential form with a temperature-dependent stretch- 
ing exponent. The Ogielski form of the nonequilibrium 
relaxation is characterized by power-law behavior in the 
intermediate time scale and faster decay in the long-time 
region. The emergence of power-law behavior in the in- 
termediate time scale suggests the presence of a sort of 
coarsening dynamics with dynamic scaling, which is sup- 
ported by the evolution of vortex patterns. 

This is also consistent with the vortex configurations 
in low-energy states obtained from a global optimiza- 
tion algorithm such as conformational space annealing 
(CSA) j 21 i 22 Those vortex configurations obtained from 
CSA exhibit interesting features: There exist approxi- 
mately parallel domains of local vortex order correspond- 
ing to the staircase ground state of / = 2/5. 23 Those 
domains are separated by domain- wall regions, consist- 
ing of characteristic local arrangements of four consecu- 
tive vacancies along one (horizontal or vertical) direction, 
thereby neighboring domains are parallel-shifted by two 
horizontal (or vertical) lattice units. Reminiscent of the 
smectic order in a liquid crystal, such configurations may 
be described as smectic (liquid-crystalline) arrangements 
of diagonal chains of vortices. Similar features of the vor- 
tex lattice were also reported for the case of frustration 
/ = 13/34 and 21/55, which are rational approximants 
to the irrational value / = 1 — 

Evolution of the vortex pattern in general exhibits 
growth of local vortex order with time toward the low- 
temperature anisotropic state described above. This 
power-law behavior of coarsening is expected to be in- 
terrupted by long-time fluctuations of domain walls, lim- 
iting the maximum local domain size to around 12 lattice 
units on average. In the long-time limit, vortex motions 
occur predominantly in the domain-wall region through 
transfer of vortices between neighboring diagonal vortex 
chains, which correspond also to the motion of four con- 
secutive vacancies. Accordingly, the higher mobility of 
the vortices located in domain- wall regions can naturally 
explain dynamic heterogeneity in the system, leading also 



to the stretched exponential relaxation. 

As for the resistive transition, one can argue that, as 
long as the domain-wall defects keep fluctuating with un- 
bounded displacement, the system would remain resistive 
(i.e., exhibiting finite resistance) even at low tempera- 
tures below the transition to the anisotropic phase, since 
there always exists a direction along which the helicity 
modulus vanishes. Those domain- wall defects are likely 
to freeze at a much lower temperature which vanishes 
in the limit of irrational frustration; thus concluded is a 
zero-temperature resistive transition in the system with 
irrational frustration. Note, however, that this argument 
applies to a pure system with no quenched disorder. In 
a real Josephson-junction array, disorder is inevitable in 
the distribution of critical currents of individual junc- 
tions, which leads to pinning of domain-wall defects at 
finite temperatures. This may explain the recent ex- 
perimental results reporting a finite-temperature resis- 
tive transition near incommensurability of the magnetic 
frustration^ 

This paper is organized as follows: In Sec. II we intro- 
duce the equations of motion for the RS J dynamics of the 
system in the fluctuating twist boundary conditions. Sec- 
tion III presents the results of simulations performed on 
the equations of motion. Both equilibrium and nonequi- 
librium relaxation behaviors of chirality autocorrelations 
are examined and coarsening, interrupted by fluctuating 
domain-wall defects, is addressed. Finally, a summary is 
given in Sec. IV. 



II. EQUATIONS OF MOTION 

We begin with the set of equations of motion for the 
phases {4>i\ of the superconducting order parameters in 
an L x L square array. In the RSJ model under the 
fluctuating twist boundary conditions, 24 they read: 



dt 



sm 



- r« • A) + C k 



= 0, 



(1) 



where we have employed the abbreviations cftij = 
<j)i—<j)j~Aij and ry = rj— Vj, and the primed summation 
runs over the nearest neighbors of grain i. The position 
of grain i is represented by = (xi,yi) with the lattice 
constant set equal to unity while the gauge field Ay is 
given by the line integral of the vector potential A: 



4>o 



A-rfl 



(2) 



with the flux quantum $ = hc/2e. The frustration pa- 
rameter /, which measures the number of flux quanta per 
plaquette, is given by the directional sum of the gauge 
field around a plaquette: 



/ 



(3) 
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In Eq. ([T]) the energy and the time have been expressed 
in units of HI c /2e and h/2eRI c , respectively, with single- 
junction critical current I c and shunt resistance R. The 
thermal noise current Qj is assumed to be white, satisfy- 



ing 



{Cij(t+rKki{t))=2k B TS(T)(S ik 6 jl -Su6 j k) (4) 



at temperature T. Henceforth we set the Boltzmann 
constant ks = 1, thus measuring the temperature in 
units of HI c /2ekB- The dynamics of the twist variables 
A = (A x ,A y ) is governed by 



dA a 1 ^ . a 

(«>« 



dt 



A a )+Ca, 



(5) 



where ^2uj\ denotes the summation over all nearest- 

x,y) and ( a 



neighboring pairs in the a-direction (a 
satisfies 



(Ca(t+T)Cb(t)) 



2T 
U 



SabS{r). 



(6) 



To study the relaxation of the system, we let the 
system evolve from some random initial configurations 
quenched to given temperatures and measure the chiral- 
ity autocorrelation function: 



Cq (t J rt w , t w ) 



with the waiting time t w . Here the chirality is defined to 
be 



Tij • A(t) 



(8) 



where X]_p denotes the directional plaquette summation 
of links around dual lattice site R and the phase differ- 
ence 4>ij(t) — ■ A(£) is defined modulo 2ir in the range 

(-7T,7r]. 

In numerical simulations, we have integrated directly 
the equations of motion (|TJ) and ([5]) via the modified Euler 
method with time step At = 0.05. The time step has 
been varied, only to give no appreciable difference. We 
have considered mostly systems of linear size L = 55 and 
89, taking averages typically over 100 to 600 ensembles 
with random initial states. The size has been chosen as 
members of the Fibonacci sequence, thus to minimize the 
boundary effects due to the irrational frustration / — 
1 — g. For comparison, we have also considered close 
rational approximants / = 21/55 and 34/89, to find no 
qualitative difference in relaxation dynamics. 



III. SIMULATION RESULTS 

We first consider the equilibrium relaxation behavior 
of the system, obtained in the following way: Dynamic 



simulations are performed with random initial states and 
then the autocorrelation function in Eq. |(7J) is computed 
for different values of the waiting time. When the wait- 
ing time is sufficiently large (and the temperature is not 
too low) , the autocorrelation function no longer depends 
on the waiting time, collapsing onto a single relaxation 
function. This collapsed relaxation function is taken to 
be the equilibrium relaxation function C eq (t). In this 
way, equilibration of the system is achieved here, down 
to temperature T = 0.14. In Fig. [IJa), the behavior 
of C eq (t) is exhibited for several values of the tempera- 
ture T. Excluding the earliest time regime, one may fit 
this relaxation behavior to a stretched exponential form: 
C eq (t) w Aexp[-(i/r)' 3 ]. Figure QJb) shows that the 
stretching exponent f3 decreases as the temperature is 
lowered, reaching the value (3 w 0.5 at T = 0.14. On 
the other hand, Fig.fljc), plotting the relaxation time r 
versus the inverse temperature T _1 , discloses the Vogel- 
Tamman-Fulcher behavior: 



t(T) = r exp 



D- 



T 



(9) 



with the fragility parameter D = 10.2 and other parame- 
ters To = 3.65 and To = 0.06. (Here To is merely a fitting 
parameter, perhaps not associated with a transition.) 
These results are consistent with those from Langevin 
dynamics 

As the temperature is further lowered, especially below 
0.14, the system, starting from a random initial state, 
does not relax to the equilibrium within the available 
computing time. Instead of pursuing equilibrium relax- 
ation at these low temperatures, we probe nonequilib- 
rium relaxation for various waiting times by letting the 
system evolve from random initial states. It is found 
that the resulting relaxation of chirality autocorrelations 
proceeds in three stages. Figure[2ja) shows the time evo- 
lution of C q {t + t w ,t w ) for the waiting time t w — 30, 000 
at temperature T = 0.15,0.13, and 0.10. For the same 
data, plotted in Fig.[2Jb) is 1 — C q (t + t w , t w ) versus time 
t, where one can see that, for almost three decades be- 
ginning from the earliest time, the relaxation fits nicely 
to the von Schweidler form 1 — (i/ro) b with b w 1/2. 
In Fig. [2jc) we show 1 — C q (t + t w ,t w ) versus time 
t at temperature T — 0.13 for various waiting times. 
It is observed that the von Schweidler behavior with 
b w 1/2 is rather robust, hardly depending on the wait- 
ing time or the target temperature. This is in contrast 
with the results of Langevin simulations, where the value 
of b tends to deviate from 1/2, getting smaller at low 
temperatures^ We presume that this discrepancy arises 
from the difference in vortex dynamics (and diffusion) at 
short times between RSJ dynamics and simple Langevin 
dynamics. 

We now turn to the relaxation behavior at intermedi- 
ate and late stages for various waiting times and tem- 
peratures. Figure [3] shows the nonequilibrium relaxation 
of the chirality autocorrelation function for null waiting 
time at temperature T = 0.15,0.14,0.13, and 0.12. It is 
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observed that the relaxation exhibits power-law behavior 
in the intermediate-time regime, followed by faster decay 
in the long-time regime. Here it is tempting to fit the 
relaxation of the autocorrelation function to the Ogielski 
form: 

C^^A^expl-it/^f] (10) 

with the exponents a and (3 given in Fig. |4]and Table I. 
Note that a depends substantially on the temperature, 
varying in the range of 0.19 to 0.35 at temperatures be- 
tween 0.08 and 0.15. 



T 


a 


P 


0.08 


0.19(1) 




0.09 


0.23(1) 




0.10 


0.258(15) 


1.27(15) 


0.11 


0.264(10) 


0.80(3) 


0.12 


0.300 (15) 


0.67(3) 


0.13 


0.300 (13) 


0.54(2) 


0.14 


0.350(10) 


0.83(4) 


0.15 


0.350(10) 


0.83(4) 



TABLE I: Exponents a and /3 depending on the temperature 
T. Note that (3 is not shown for the cases of T = 0.08 and 0.09. 
At these low temperatures, relaxation is too slow to observe 
clearly the late-time stretched exponential part within the 
computational time window, making it formidable to estimate 
f3. 

Also shown in Fig [5] is the nonequilibrium relaxation 
of the chirality autocorrelation function for various wait- 
ing times, at temperature T = 0.14,0.13, and 0.12. 
For longer waiting times [t w > 100.0), this behavior 
crosses over to the stretched exponential form C q (t) sa 
A2 exp [-(t/r 2 ) /3 ]. Emer gence of the power-law behavior 
of the nonequilibrium relaxation at intermediate times 
strongly suggests that there exist some coarsening pro- 
cesses in the system. Figure [5] exhibits snapshots of the 
vortex configuration at temperature T = 0.13, taken at 
several time instants; one can recognize slow growth of lo- 
cal order corresponding to / = 2/5 vortex patterns, with 
diagonal chain structures. Since the frustration of the 
system is given by / = 1— g, which is slightly smaller than 
2/5, these locally ordered domains of / = 2/5 patterns 
may not grow to span the whole system. Instead, there 
should exist finite length (and also time) scales for the 
growth of these local domains, beyond which the growth 
is interrupted by domain- wall regions of lower vortex den- 
sity, so that the net vortex density of the whole system 
becomes equal to 1 — g. 

In order to understand the vortex configuration at- 
tained in the long-time limit at low temperatures, we in- 
vestigate the configuration of low-energy states by means 
of the efficient optimization algorithm, CSAi^l Snapshots 
of typical vortex configurations obtained via CSA are 
shown in Fig where we observe domains of locally 
ordered vortex patterns corresponding to the staircase 



ground state of the system with frustration (or vortex 
density) / = 2/5. These domains of typical width 9 to 
12 lattice spacings are separated by domain-wall (line 
defect) regions that consist of characteristic local ar- 
rangements of four neighboring vacancies (see bar-shaped 
regions, each with four consecutive empty plaquettes); 
there the vortex density is lower than the locally ordered 
regions (of local vortex density 2/5) in such a way that 
the net vortex density of the whole system is precisely 
equal to 1 — g (which is less than 2/5). This config- 
urations may also be described as a liquid-crystal-type 
arrangement of diagonal chains of vortices (of length 9 
to 12), with the neighboring chains of vortices shifted in 
the diagonal direction by about half the length of the 
chains. Similar configurations were observed in the case 
of rational approximants to 1 — g 

In the long-time limit, we expect that the coarsening 
dynamics will lead eventually to the locally ordered con- 
figuration with domain walls, shown in Fig. [7J Within 
a locally ordered domain, vortices are almost rigid and 
resist moving. In contrast, those vortices at ends of the 
chains are easily put in motion, jumping into the bar- 
shaped vacancies in domain-wall regions and thus join- 
ing another vortex chain. This in turn gives rise to 
domain- wall fluctuations (or equivalently, fluctuations in 
the length of the vortex chains), which would lead slowly 
but ultimately to the complete restructuring of the local 
vortex configuration. It is thus expected that dynamic 
heterogeneity naturally emerges from the existence of dis- 
tinctly mobile vortices in the domain- wall regions 

We believe that such domain-wall fluctuations and re- 
structuring of the local vortex configuration can explain 
the absence of freezing of the relaxation at low tempera- 
tures. Namely, the low-temperature state is presumably 
of a liquid crystalline type: Even though there exists ori- 
entational order in the chain-like arrangement of vortices, 
those diagonal chains can flow (like a liquid) due to the 
transfer of vortices between neighboring chains, corre- 
sponding to the fluctuating motion of the bar-shaped va- 
cancies and resulting in finite resistivity. At a much lower 
temperature one may expect complete freezing of the do- 
main walls to occur. The freezing temperature should 
depend on the commensurability of the distribution of 
the vacancy defects and the underlying background vor- 
tex lattice; this is directly related to the rationality of 
f = p/q (with p and q relatively prime integers) and the 
freezing temperature is expected to vanish in the irra- 
tional limit (q — > 00). 

It is of some interest to note the similarity to the 
behavior found in the lattice coulomb gas with charge 
density near the golden number^ Even though the de- 
tailed ordering pattern is different, it was found numer- 
ically that there exist two-step transitions in the lattice 
coulomb gas with 1/3 < / < 2/5, where the intermediate 
phase corresponds to anisotropic striped charge order- 
ing together with mobile charges within partially-filled 
channels. There the lower transition corresponds to the 
complete freezing of charges within partially filled diag- 
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onal channels, which occurs at temperatures sensitively 
dependent on the rationality of / due to the commensu- 
rability effects. In this study our system is the Josephson- 
junction array in the limit of irrational frustration, and 
we thus expect that the true vortex freezing would occur 
at zero temperature. 

As long as the domain-wall defects fluctuate with un- 
bounded displacement, the system should remain resis- 
tive, exhibiting finite resistance. These domain-wall de- 
fects are expected to freeze at much lower temperature 
of the order 1/q for / = p/q; this leads to the zero- 
temperature resistive transition in the limit of irrational 
frustration. Note that this argument applies for pure 
systems with no quenched disorder. In the case of real 
Josephson-junction arrays, some type of disorder, e.g., in 
the distribution of the critical currents of individual junc- 
tions is unavoidable. Therefore pinning of the domain- 
wall defects can easily occur at finite temperatures, hav- 
ing resistance vanish. This may explain the recent ex- 
periment reporting a finite-temperature resistive transi- 
tion near incommensurability of the magnetic frustration, 
where disorder in the critical currents of individual junc- 
tions was noted to be up to 15%^ 

IV. SUMMARY 

We have studied the relaxation behavior of the chi- 
rality autocorrelation function in two-dimensional super- 
conducting arrays under irrational frustration at inter- 
mediate and low temperatures. Both equilibrium relax- 
ation and nonequilibrium relaxation have been investi- 
gated via numerical simulations of RSJ dynamics. Equi- 
librium relaxation dynamics reveals characteristic fea- 



tures of the stretched exponential form with the Vogel- 
Tamman-Fulcher dependence of the relaxation time. 

Nonequilibrium relaxation at low temperatures, be- 
ginning with random initial states exhibits interest- 
ing waiting-time dependence: For short waiting times, 
the relaxation follows a power-law behavior in the 
intermediate-time regime and faster decay in the long- 
time regime, characterized together by the Ogielski form 
with temperature dependent exponents. For longer wait- 
ing times, this gradually crosses over to a stretched expo- 
nential form. Further, in all cases of short and long wait- 
ing times, the relaxation at early time stage fits nicely to 
the von Schweidler form with exponent about 1/2. 

It has been argued that the power-law behavior orig- 
inates from coarsening dynamics up to a certain length 
scale, with the local vortex order corresponding to the 
ground state of / = 2/5, found to be consistent with the 
vortex configurations of low-energy states. This coarsen- 
ing, however, is interrupted by the presence of domain- 
wall defects. Such chain-like domain- wall regions are ex- 
pected to provide possible mechanism of dynamic het- 
erogeneity and stretched exponential relaxation. Since 
the RSJ dynamics adopted in this study can be realized 
experimentally, it would be of interest to carry out ex- 
periment and compare the results. 
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FIG. 1: (a) Relaxation of the chirality autocorrelation func- 
tion (at equilibrium) C eq (t) at various temperatures T in the 
irrationally frustrated array of linear size L = 89. Stretched 
exponential fits are also plotted, (b) Stretching exponent ver- 
sus the temperature. The dotted line is merely a guide to the 
eye. (c) Relaxation time versus the inverse temperature. The 
dashed line represents the Vogel-Tammann-Fulcher fit with 
To = 0.06 and to = 3.65. The inset shows the result in a 
different scale. 
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FIG. 2: (a) Relaxation of the chirality autocorrelation func- 
tion C q (t+t w , t w ) with time t in the early-time regime at tem- 
peratures T = 0.10,0.13, and 0.15 for L = 89. The waiting 
time is chosen to be t w = 30000. (b) The log-log plot of 
1 — C q (t+t w ,t w ) versus t for the same data. For almost three 
decades in the early-time regime, the relaxation behavior is 
observed to follow 1 — C q (t+t w , t w ) ~ t 1 ' 2 at all the three tem- 
peratures, as indicated by the dot-dashed line, (c) The log-log 
plot of 1 — C q (t+t w ,t w ) versus t at temperature T = 0.13 for 
three different waiting times t w = 1.0,100 and 30000. It is 
shown that the early-time behavior does not depend qualita- 
tively on the waiting time. 
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FIG. 3: Nonequilibrium relaxation of the chirality autocor- 
relation function C q (t+t w , t w ) for zero waiting time (t w = 
0.0), starting from random initial states, at temperature 
T = 0.15, 0.14, 0.13 and 0.12. While in the intermediate-time 
regime the nonequilibrium relaxation function exhibits fea- 
tures of power-law behavior, this behavior is interrupted by 
faster decay in the long-time regime, which together can be 
fitted to the Ogielski form, as shown by the dot-dashed lines. 
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FIG. 4: Exponents (a) a and (b) (3 versus temperature T in 
the fit of the non-equilibrium relaxation for t w = to the 
Ogielski form. 
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FIG. 5: Aging relaxation of the chirality autocorrelation func- 
tion C q (t+t m ,t m ) for various values of the waiting time t w = 
0, 1, 3, 10, 30, 100, 300, 1000, 3000, 10000, and 30000 (from left 
to right), starting from random initial states, at temperature 
T = (a) 0.14, (b) 0.13, and (c) 0.12. Power- law behavior t~ a 
in the intermediate-time regime is manifested by dotted lines 
with the exponent a = (a) 0.35, (b) 0.33, and (c) 0.31. As the 
waiting time grows, the relaxation develops simple stretched 
exponential behavior. 
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FIG. 6: Snapshots of the vortex configuration in the system 
of size L — 55, beginning from a random initial state and 
quenched to temperature T = 0.130, at time t = (a) 0.0, (b) 
8.0, (c) 128, and (d) 2048. 
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FIG. 7: Snapshots of the vortex configuration obtained via 
CSA for lattice size L = 55 with energy E = (a) -1.273786 
and (b) -1.274077 (in units of hl c /2e), respectively. 



